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ABSTRACT 

We present the results of an analysis of a well-selected sample of galaxies with 
active and inactive galactic nuclei from the Sloan Digital Sky Survey, in the range 
0.01 < z < 0.16. The SDSS galaxy catalogue was split into two classes of active galax- 
ies, Type 2 AGN and composites, and one set of inactive, star- forming/passive galax- 
ies. For each active galaxy, two inactive control galaxies were selected by matching 
redshift, absolute magnitude, inclination, and radius. The sample of inactive galaxies 
naturally divides into a red and a blue sequence, while the vast majority of AGN 
hosts occur along the red sequence. In terms of Ha equivalent width, the popula- 
tion of composite galaxies peaks in the valley between the two modes, suggesting a 
transition population. However, this effect is not observed in other properties such 
as colour-magnitude space, or colour-concentration plane. Active galaxies are seen to 
be generally bulge-dominated systems, but with enhanced Ha emission compared to 
inactive red-sequence galaxies. AGN and composites also occur in less dense environ- 
ments than inactive red-sequence galaxies, implying that the fuelling of AGN is more 
restricted in high-density environments. These results are therefore inconsistent with 
theories in which AGN host galaxies are a 'transition' population. We also introduce 
a systematic 3D spectroscopic imaging survey, to quantify and compare the gaseous 
and stellar kinematics of a well-selected, distance-limited sample of up to 20 nearby 
Seyfert galaxies, and 20 inactive control galaxies with well-matched optical proper- 
ties. The survey aims to search for dynamical triggers of nuclear activity and address 
outstanding controversies in optical/IR imaging surveys. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) have long been considered 
a curiosity in their own righ t (ISchmidtlll963l : iLvnden-Benl 
ll969l : lBlandford fc Reeslll974 ). but are now recognised to be 
integral to galaxy formation and evolution. At early cosmo- 
logical epochs, gas- rich g alaxies are thought to form and col- 
lide in violent mergers (|Mihos fc Hernquistlll996h . trigger- 
ing vast bursts of star form ation, and fuelling supermassive 
black holes in their cores (|Kriek et al.ll2006lh Recent sim- 
ulatio ns of isolated and merging galaxies by ISpringel et all 
have incorporated feedback from star formation and 
black-hole accretion, and find that once an accreting super- 
massive black hole (SMBH) has grown to some critical size, 
the AGN feedback terminates its growth as a large fraction 
of the remaining nuclear gas is driven out by the powerful 
quasar. 

In the current epoch, the peak of the quasar era is 



E-mail:pbw@astro. livjm.ac.uk 



over and the galaxy merger rate has declined (|Strucklll997h . 
Nevertheless, at least 20% of today's galaxies show scaled- 
down quasar activity in their centres, and direct measure- 
ments of active and inactive galaxy dynamics have revealed 
a tight correlation between the central black-hole mass and 
host galaxy stell ar velocity dispersion or bulge mass. This 
Mm- a relation jGebhardt et al.ll200t)l : iMerritt fc Ferraresel 
l200ll ) points to an intimate link between host galaxy evolu- 
tion and central black-hole growth, suggesting that all bulge- 
dominated galaxies today harbour dead quasars and that 
a lack of nuclear activity cannot be attributed to the ab- 
sence of a central black hole. Therefore, given the ubiquity 
of supermassive black holes, what determines the degree of 
nuclear activity in todays galaxies and what role is played 
by the host galaxy in triggering and fuelling their dormant 
black holes remain open issues. 

Previously, studies of nearby active galaxies were based 
on small galaxy sam ples, like the C/A Seyfert Sample 
dHuchra fc Burd|l992h . the 12 Micron Active Galaxy Sam- 
ple (|Rush et alJll993h and the lHo et ail (|l995l ) galaxy sam- 
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pie of approximately 486 galaxies covering a range of ac- 
tivity types. Now, however, the standard has been set by 
the Sloan Digital Sky Survey (SDSS), from which a range 
of useful samples of both active and inactive galaxies can be 
selected in a we ll-defined, unifor m way. 

The SDSS (jYork et all2000h provides, for the first time, 
5-band photometry an d spectroscopy of many thousands 
of low redshift AGN l|Kauffmann et all 12003d ; lHao et all 
l2005al ). enabling constraints to be placed on galaxy and 
AGN evolution over a wide range of galaxy masses, and 
acts as the definitive supporting data to all detailed follow- 
up galaxy studies. Initial investigations into the growth and 
evolution of black holes using these databases have yielded 
contrasting results — the origin of the correlation between 
galaxy bulge and central black-hole masses is hotly debated 
in p articular. 

iHeckman et al.l l|2004h find that the majority of present- 
day accretion occurs onto 10 s solar-mass black holes in 
moderate mass galaxies, suggesting that bulge and black- 
hole evolution is still tightly coupled today, and that the 
evolution of AGN luminosity functions is driven by a de- 
crease in the mass s c ale of accreting black holes. In con- 
trast, the lHao et "all (l2005bl ) study of about 3 000 SDSS 
AGN concludes that evolution in AGN luminosity func- 
tions is driven by evolution in th e Eddington ratio, rather 
than black-hole mass. Mean while, iGrupe fc Mathuri (|2004l ) 
and iMathur fc Grupd (|2005h use a sample of 75 X-ray se- 
lected AGN to argue that black holes in Narrow Line Seyfert 
1 galaxies in particular, grow by accretion in well- formed 
bulges to produce the M. — a relation over time, refut- 
ing theories for the origin of black-hole / bulge relations in 
which black-hole mass is a constant fraction of bulge mass 
at all epo chs or in which bulge growth is controlled by AGN 
feedback |King||2003l ; rFerraresdl2002l ; lKriek et aLll2006l ). The 
su bclass of NLSls has be en further explored more recently 
bv lKomossa fc Xul {2007), who find that NLSls are accret- 
ing at a rate higher than the Eddington rate, confirming that 
their BHs must be growing. They suggest that either NLS1 
galaxies evolve into Broad Line Seyfert 1 (BLS1) galaxies 
with respect to their black hole mass distribution, which 
would require some change in the bulge properties, possibly 
due to feedback, or that NLSls are just low- mass extensions 
of BLS1 galaxies, and the high accretion rate could just be 
caused by a relatively short-lived accretion phase. 

Another characteristic revealed by large automated sur- 
veys such as the SDSS, is t hat the galaxy population is found 
to be bimodal in colour |Lillv et alj [l995l ; IStrateva et al.l 
2001), with it now being more natural to describe a galaxy 
as being on the "red sequence" or "bl ue sequence", rathe r 
than being "early type" or "late type" (|Baldrv et al.| [2006). 
A key goal of galaxy evolution theory is then to explain the 
colour bimodality of galaxies, the relationships within each 
sequence, and where active hosts fit into this picture. There- 
fore, the understanding of galaxy formation and evolution 
processes necessitates full inclusion of AGN and their hosts. 

Recently efforts have been made to try to understand 
where AGN host galaxies fall in colour-magnitude space. 
The red sequence consists of mainly massive, passively 
evolving galaxies, while the majority of galaxies show blue 
colours, attributed to ongoing star formation. The two se- 
quences are separated by a relatively narrow valley in colour 
space. 



The emerging consensus is that intense star formation 
is fuelled by galaxy mergers at high redshift, which forms 
massive bulges, and at some point the star formation ceases, 
resulting in the galaxy migrating from the blue sequence to 
the red sequence. AGN have been singled out as the mech- 
anism for star-formation quenching, leading to suggestions 
that AGN should occupy a distinct, 'transition' region, of the 
colour- magnitude diagram (CMD). However this does not 
necessarily have to be the case, as the resulting galaxy could 
just be bluer due to increased star-burst activity, or redder 
due to enhanced dust (e.g., Luminous Infra- Red Galaxies), 
or some combination of the original colours of the merging 
galaxies. Young stellar populations and nuclear starbursts 
are therefore also important components in the dynamics 
of nu clear activty l|Gonzalez Deleado et al.l l 1998 : Sarz i et al.l 
120071) . 

INandra et al l (|2007t ) studied the colour-magnitude rela- 
tion for a sample of 50 X-ray selected AGN from the AEGIS 
survey (All-wavelength Extended Groth strip International 
Survey), in the range 0.6 < z < 1.4. They conclude that 
AGN fall on the red sequence or on the red edge of the blue 
sequence, with many in between these two modes. 

iMartin et alj (|2007f ) further explored the idea of a 'tran- 
sition' region by exploiting the separation of the blue and red 
sequences in the UV — r colour-magnitude diagram. Using 
a sample of UV selected galaxies from the GALEX Medium 
Imaging Survey, along with SDSS data, they explored the 
nature of galaxies in the transition zone. The AGN fraction 
of their sample was found to peak in the transition zone, and 
they also found circumstantial evidence that star-formation 
quenching rates were higher in higher luminosity AGN. 

Higher image quality than the SDS S is provided by 
the Millennium G alaxy Catalogue (MGC; iLiske et al.1l2003l ; 
ICross et alll2004l ). which has obtained redshifts for 10 095 
galaxies to B < 20 mag, covering 37 deg 2 of equatorial sky. 
The MGC study of galaxy bimoda lity was in the colour- 
structure plane, (I Driver et al.ll2006l ; hereafter D06) and ex- 
pressed a contrasting theory for the bimodality. The MGC 
survey has suffi cient resolution to achieve reliable bulge-disc 
decomposition ijAllen et al]l2006l ). and in separating out the 
two components, D06 claim that galaxy bimodality is caused 
by the bulge-disc nature of galaxies, and not by two dis- 
tinct galaxy classes at different evolutionary stages. In this 
case, they show that the bulge-dominated, early-type galax- 
ies populate one peak and the bulge-less, late-type galaxies 
occupy the second. The early- and mid-type spirals sprawl 
across and between the peaks. They also propose that the 
reason for the dual structure of galaxies is due to galaxy 
formation proceeding in two stages; first, there is an initial 
collapse phase, which forms the centrally concentrated core 
and a black hole, and second, there is the formation of a 
planar rotating disc caused by accretion of external mate- 
rial building up the galaxy disk. 

In this paper, we present the properties of a magnitude- 
limited sample of "active" and "inactive" galaxies carefully 
selected from the SDSS, which also acts as the parent sample 
for detailed followup studies using the IMACS-IFU on the 
Magellan 6.5m telescope. Our sample selection is described 
in §3, the sample properties in §4, and implications of these 
properties in §5. We conclude by introducing the Magellan 
survey, which is now underway. 
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2 THE SLOAN DIGITAL SKY SURVEY 

The SDSS |York et all l200d : IStoughton etafl 120021 1 is an 

imaging and spectroscopic survey, which, when completed, 
will provide detailed optical images covering more than a 
quarter of the sky and a 3-dimensional map of about a 
million galaxies and quasars. The survey uses a dedicated 
2.5m telescope at Apache Point Observatory, New Mexico, 
equipped with a large format mosaic CCD camera — 30 
CCDs in six columns and five rows — to image in five op- 
tical bands; u, g, r, i and z, with eff ective wavelengths of 



Solid contours: Nil 8/N > 2-0 



Dashed contours: 



1 19961 ; 



3 550, 4 670 6160 , 7 470 and 8 920 A (|Fukugita et al 
iGunn et al.ll 19981 ). 

The images are calibr ated photometri cally l|Hogg et al.l 

l200ll ) and astrometrically dPier et al.ll2003l ) using a separate 
0.5-m telescope, and reduced using a pipeline photo that 
measures the observing conditions, and detects and mea- 
sures objects. In particular, photo produces various types 
of magnitude including 'Petrosian', the summed flux in 
an aperture that depends on the surface-brightness pro- 
file of the object. The magnitud es are Galactic extinction- 
corrected using the dust maps of lSchlegel et all dl99Sl). De- 
tails o f the imaging pipelines are given by Istoughton et al.l 
|2002| y 

The main telescope is also equipped with two double 
fibre- fed spectrographs, covering a wavelength range of 3800- 
9200 A, over 4 098 pixels, with a resolution of A/ A A ~ 1 800. 
Once a sufficiently large area of sky has been imaged, the 
data are analysed using 'targeting' software routines that 
determine the objects to be observed spectroscopically. The 
targets are then assigned to plates, each with 640 fibres, 
using a tiling algorithm l|Blanton et al.l l2003al ). The main 
restriction is that two fibres cannot be placed within 55" on 
the same plate. 

Spectra are typically taken using three 15- minute expo- 
sures in moderate conditions (the best conditions are used 
for imaging). The signal-to-noise ratio (S/N) is typically 10 
per pixel (1-2A) for galaxies. The pipeline spec2d extracts, 
and flux- and wavelength-calibrates the spectra. 



3 SAMPLE SELECTION 

We start with the gal axy catalogue from the SDSS Fo urth 
Data Release (DR4; lAdelman-McCarthv et all 120061 ). as 
distribu ted by the MPA Gar c hing group. The cata- 
logued d Kauffmann et all l2003al ; iBrinchmann et aL I I2004 
iTremonti et all |2004| ) contains photometric and follow-up 
spectroscopic data for all objects spectroscopically classi- 
fied as galaxies with Petrosian r-band magnitudes less than 
17.77. The catalogue does not include objects classified 
as 'quasars' by the SDSS spectral classification algorithm. 
Broad emission-line galaxies ('Type 1' Seyferts) are there- 
fore generally excluded from the MPA catalogue, although 
there are still some Type 1 Seyferts remaining if their spec- 
tra contain a significant fraction of stellar light. Galaxies 
with a velocity dispersion, a, of greater than 500 km s _1 are 
therefore assumed to be broadline Seyferts (Type 1 AGN), 
and due to the relatively low number of these (less than 4 




Figure 1. BPT diagnostic diagram for DR4 SDSS galaxies, show- 
ing the distributions of the sample, above and below our S/N cut 
of 2.0. The dot-dashed line is the KewOl maximum star-burst 
line. The dotted line is the Kauf03 classification line. 



per cent of the overall sample), are rejected from this anal- 
ysis. In this paper therefore, we will only be comparing in- 
active galaxies with narrow-line galaxies ('Type 2' AGN, or 
'Type 2' Seyferts). Exclusion of type 1 AGN from our active 
sample does not affect the overall conclusions. 

The MPA DR4 catalogue contains derived physical 
properties for 567 486 galaxy spectra of 520 738 individual 
galaxies, with duplicates being removed using a 1.5" match- 
ing radius. A redshift cut of 0.01 < z < 0.16 is then ap- 
plied to all galaxies to limit evolution, k-corrections, and 
low-redshift problems. 359 154 DR4 galaxies then remain in 
our sample, from which active galaxies are isolated, leaving 
most of the remaining galaxies to be used as potential con- 
trol galaxies. Once the AGN sample has been produced, two 
inactive control galaxies are selected for each AGN, closely 
matched in their optical properties. 



3.1 AGN and Inactive Classification 

In order to compare active and inactive galaxies, the SDSS 
galaxies must be cl assified and the AGN separated out. 
iBaldwin et ail l|l98ll) . hereafter BPT, demonstrated that it 
is possible to distinguish Type 2 AGN from normal star- 
forming galaxies using the relative intensity ratios of easily- 
measured emission lines. BPT found that the most use- 
ful intensity ratios were [OIII]/IL3, [Nil] /Ha, [SII]/Ha and 
[OI]/Hq, and plotting these ratios against each other yields 
diagnostic diagrams. 

These emission lines have be en accurately measure d for 
all DR4 galaxies, as described by ITremonti et al.l l2004l. The 
SDSS spectroscopic pipeline performs only a simple estimate 
of the stellar continuum, which is good for strong lines, but 
not for weaker emission lines. The Tremonti et al. pipeline 



on the other hand, estimates the stellar continuum by taking 
1 Data publicly available at http://www.mpa-garching.mpg.de/SDSS/DMo' account stellar Balmer absorption, and hence recovers 
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weaker nebular emission line features. This allows for a more 
accurate classification scheme of SDSS galaxies. 

We adopt a similar clas sification scheme to that devised 
bv lBrinchmann et al.l (|2004T l . who split the sample into three 
main classes of galaxies based on their positions in the var- 
ious BPT diagrams. These classes were defined to be star- 
forming or passive galaxies, AGN, and composite galaxies. 
Two extra sub-samples were derived for low signal-to-noise 
star-forming galaxies, and low signal-to-noise AGN. A final 
class was designated to ambiguous galaxies whose position 
in the various BPT diagrams appeared to place them in 
more then one class - this class most likely to be made up 
of galaxies with weak or no emission lines. 

For our purpose, it is sufficient to separate 'active' 
and 'inactive' galaxies. We therefore split the sample into 
three classes - star-forming/passive gala xies, Type 2 AGN, 
and composite galaxies. We also use the lKewlev et af1l200ll 
(hereafter KewOl) maximum star-burst separation line (dot- 
dashed line in Fig.[TJ and the Kauffmann et al. 2003bl (here- 
after Kauf03) demarcation line (dotted line in Fig. [1} to se- 
lect our sample. Therefore, in this study, a galaxy is defined 
to be a Type 2 AGN if 

\og([OIII]/HB) > 0.61/(log([A7/]/#a) - 0.05) + 1.3 (1) 

or defined to be a purely star-forming galaxy if 

\og([OIII]/HB) < 0M/Qog([NII]/Ha) - 0.47) + 1.19. (2) 

where [OIII] refers to the [OIII] A5007 emission line, and 
[Nil] is the [Nil] A6584 lineQ Galaxies that fall in between 
these classification lines are galaxies that are classified as 
AGN according to the more lenient Kauf03 separation line. 
These galaxies are kept as the separate class of composite 
galaxies, as per Kauf03. 

We tested the sensitivity of our results to the galaxy 
classification, by slightly altering equations (1) and (2) in 
our selection process, making the composite sample smaller. 
The effects of such changes were minimal, so for simplicity, 
and comparison with other results, we continued to use the 
KewOl and Kauf03 separation lines. 

The Type 2 AGN and composite galaxies were also sub- 
jected to a signal-to-noise cut, to ensure reliable AGN clas- 
sification. We chose to cut at S/N > 2.0 in the [Nil] A6584 
emission line. As the [Nil] emisson line flux for star-forming 
and passive galaxies can be low, no S/N cut was required for 
that classification. Figure [1] shows a BPT plot for the over- 
all DR4 sample. We have plotted the sample in two S/N 
bins. The solid contours show the distribution of galaxies 
with S/N > 2.0, while the dotted contours show the dis- 
tribution of those with S/N < 2.0. The latter distribution 
shows that these low S/N galaxies are randomly distributed 
across the whole BPT diagram, so will not accurately repre- 
sent AGN. Above a S/N of two, the BPT diagram takes its 
familiar form, so although this is a loose constraint it is still 
a reasonable cut for our statistical purposes in this paper, 
and allows us to keep the weaker AGN. A stricter cut of 
S/N > 3.0 was also tried, and found to have little impact 

2 When calculating the logarithms in Eqns. [l]and[2] for each line 
flux, the value used is the higher of the measured line flux and 
its 1-sigma uncertainty. This avoids anomalously large, small or 
negative line ratios when one or both lines are undetected. 



Table 1. Basic sample data following classification 



Subsample 


Number 


Percentage 


All galaxies (0.01 < z < 0.16) 


359 154 


100.0 


SF + Passive 


176 596 


49.2 


Type 2 AGN (S/N > 2 in [Nil]) 


69 783 


19.4 


Composites (S/N > 2 in [Nil]) 


64136 


17.9 


Type 1 AGN (reject) 


13 331 


3.7 


Unclassified (reject) 


35 308 


9.8 



on the overall results. Table \T\ summarises the numbers of 
galaxies classified and rejected for our analysis. 

Typical [OIII] luminosities in our active sample range 
from 1O 5 L to 1O 8 L . The mean luminosity is approxi- 
mately 10 6 Lq: the majorit y of active galaxies would be con- 
sidered to be 'weak AGN' (|Kauffmann et al.l 12003d ) as less 
than 10% of the active sample has L[OIII] > 10 7 Lq. 

3.2 Control Galaxy Selection 

In order to perform a comparative study of active and in- 
active host galaxies, careful selection of inactive control 
galaxies is required. Control matching is critical because the 
AGN classification is strongly dependent on aperture effects, 
which can be minimised by matching redshift, and inclina- 
tion, which can be minimised by matching on aspect ratio. 
We also want to match to galaxies of a similar size and lumi- 
nosity. Therefore, for each of the active galaxies we selected 
two inactive galaxies, closely matched in the following prop- 
erties: 

(i) Redshift, z, where 0.01 < z < 0.16 (AGN and com- 
posites restricted to 0.02 < z < 0.15); 

(ii) Absolute r-band magnitude, M r ; 

(iii) Aspect ratio, b/a, (Isophotal minor/major axis ratio 
in r-band); 

(iv) Radius, petroR90 (in arcsec), containing 90% of the 
Petrosian flux, averaged over r- and i-bands. 

Isophotal major and minor axis values, given by isoA 
and isoB respectively, are from the SDSS pipeline frames, 
frames measures the radius of a particular isophote of the 
galaxy, as a function of angle, and then Fourier expands this 
function. It then extracts isoA and isoB from the result- 
ing coefficients. This is done in all wavebands, but only the 
r-band values were used in this analysis. petroR90 is also 
obtained through the frames pipeline. This calculates the 
radius containing 90% of the Petrosian flux for each band. 
At this stage petroR50 is also derived, which is the radius 
containing 50% of the Petrosian flux. 

The r-band absolute magnitude used in this paper is 
given by 

M r = r~k r ~5 \og(D L /10 pc) (3) 

where r is the Milky- Way-extinction-corrected Petrosian 
magnitude, Dl is the luminosity distance for a cosmology 
with (n m ,fi A )o = (0.3,0.7) and H = 70 k m s" 1 Mpc" 1 
and k r is the fc-correction using the method of lBlanton et al.l 
l|2003bl ).FI 

3 The fc-corrections were derived from kcorrect v. 4.1.3. 
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Figure 2. Left-hand side: Normalised distributions of (a) redshift, (b) absolute magnitude, (c) aspect ratio, and (d) petroR90 radius, for 
our AGN (solid red lines) and composite galaxies (dashed green lines), along with their matched controls (dotted blue lines). Right-hand 
side: Normalised histograms of the individual differences between AGN/composites and each of their controls. The best controls are 
represented by the dotted lines, the 2nd best controls by the solid lines. 



Two controls were selected for an AGN by calculating 
the differences, Az, A[b/a], AM r , ApetroR90, between the 
AGN and each control. The differences were then weighted, 
and summed as follows: 



0.01 0.2 



A[b/a]\ \ApetroR90\ 



0.1 



+ 



0.2" 



(4) 



The two controls with the lowest values of AC were then se- 
lected, and added to the controls catalogue. This was done 
for all Type 2 AGN, and then done separately for all compos- 
ite galaxies. The weightings were such that we constrained 
the redshift and absolute magnitude more tightly then the 
other two parameters [see Fig. [2je-h)] . 



3.3 Sample Distributions 

The classification and control matching results in 2 pairs of 
galaxy catalogues: an AGN catalogue of 64 584 Type 2 AGN 
according to the KewOl criteria (inclusive of LINERs), along 
with a catalogue of 129168 control galaxies (Controls A); 
and a catalogue of 60 038 composite galaxies, with a matched 
set of 120 076 control galaxies (Controls B). 

The resulting overall distributions of the selection prop- 
erties are displayed in Fig.[2ja-d), along with the overall dis- 
tributions of the differences in these properties between the 
AGN/composites and each of their two controls - Fig. [2|e- 
h). The best controls (dotted lines) result in a sharper peak 
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Figure 3. Observed distribution of Ho equivalent width vs. absolute magnitude. The contours are determined for galaxy number counts 
in 0.135 Ha x 0.135 mag bins. The contour levels are on a logarithmic scale, starting at 64 and doubling each contour. The solid lines in 
(a) are Type 2 AGN, and in (b) are composites. The dashed lines are the control galaxies in both cases. 



in the distribution, while the second best controls in our se- 
lection — solid lines in Fig. 2(e-h) — show slightly broader 
wings. Control matching was done using all 'inactive' galax- 
ies each time, so it is possible the same controls can be picked 
more than once, if they fit the selection criteria for multi- 
ple AGN. It was not appropriate to make all control galaxies 
unique, as this would introduce systematic errors in the con- 
trol matching. 



4 RESULTS 

A significant feature in recent analyses of galaxy popula- 
tions is the recognition that the distribut ion of rest-frame 
colours is bimodal fe.g., lBaldrv et al]|2004 ). There is a well- 
characterised red sequence of predominantly massive, pas- 
sively evolving early- type galaxies, and a blue sequence of 
late type, star-forming galaxies. This is a key feature for un- 
derstanding galaxy evolution, and here we extend this anal- 
ysis to identify where active gal axies fit int o this p icture. 

Past observations by, e.g., iTreu et "al] (|2005l ). suggest 
that galaxies transition from the blue to the red sequence, in 
a ma nner consistent with "cosmic downsizing" l|Cowie et al.l 
1 19961 ); mergers between gas-rich blue galaxies producing el- 
lipticals, with a predicted phase of bright quasar activity 
which terminates star-formation, and hence leads to the 
galaxy becoming red. 

4.1 Equivalent width— Magnitude relations 

The Balm er Ha line is a clear in dicator of active star 
formation (|Kennicutt fc Kentj[l983 ). and hence is a good 
way of separating out the red and b lue galaxy sequences 
|Balogh et al.ll2004 ISames et al.ll2007r i. Fig.^a) shows the 
distribution of the equivalent width (EW) of the Ha lines in 
our AGN sample and the Controls A sample, as a function 
of absolute r-band magnitude. 



The galaxy bimodality is clearly observed in the dis- 
tribution of control galaxies, with late-type, star-forming 
galaxies (blue sequence) occupying the top half of the plot, 
and early-type, passive galaxies falling at low Ha strength, 
and peaking at a slightly higher luminosity. The peak of 
the AGN population falls along the blue edge of the red 
sequence, i.e., shows increased Ha EW relative to inac- 
tive red-sequence galaxies. The composite galaxies, shown 
in Fig. E^b) , however, show the peak of their distribution is 
located in the valley between the blue and red sequences of 
the Controls B sample. Both these results suggest evidence 
for a "transition" population, supporting galaxy evolution- 
ary theories, in which star formation in massive galaxies 
ceases due to the presence of an AGN, resulting in a mi- 
gration from the blue sequence to the red sequence. (Note 
that there is no requirement that AGN or composites have 
significant Ha emission, only Nil S/N > 2 is required). 

4.2 Colour-Magnitude Relations 

g — r colour reflects a longer timescale of star formation 
compared to Ha, which is more representative of instan- 
taneous star formation (or AGN-heated nebular emission). 
The colour-magnitude diagrams for our SDSS sample are 
shown in Figure [4] The left-hand side plot again shows the 
Type 2 AGN and Controls A, while the right-hand panel 
shows the composite galaxies with Controls B. The bimodal- 
ity of colours can be seen in the distributions of control 
galaxies in both cases. From the CMD in Fig. |4ja), the 
majority of AGN now appear to reside along the peak of 
the red-sequence, rather than along its blue-edge. The AGN 
population does, however, still extend to the red edge of the 
blue sequence as well. 

The location of the composite galaxies is less conclusive 
in this case [Fig. [4jb)] . The peak of this population is not so 
clearly just in between the red and blue sequences, as was 
the case in the Ha plot, but rather shows a more extended, 
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Figure 4. Observed distribution of colour vs. absolute magnitude, for (a) Type 2 AGN, and (b) composite galaxies, and their respective 
controls. The contours are determined for galaxy number counts in 0.1 colour X 0.06 mag bins. The contour levels are on a logarithmic 
scale, starting at 16 and doubling each contour. The solid lines represent the 'active' galaxies, while the dashed lines are the control 
galaxies. 



irregular distribution. The distribution tracks the red se- 
quence, in a manner similar to that of the AGN sample, 
but also extends to much bluer colours, indicating ongoing 
star formation, but with AGN activity also taking place. 
This result suggests that the composites do not necessarily 
represent a transition population, as implied by Fig. O^b). 

4.3 Colour-Concentration Relations as a function 

Of M r 

Although we have demonstrated that the galaxy population 
separates into red and blue sequences, this does not take into 
consideration the structure or morphology of galaxies. A ro- 
bust method of outlining the seque nces is to consider a joint 
distribution in colour and structure (I Driver et al .2006). Fig- 
ure [S] shows the distribution of observed galaxies in colour 
versus concentration index for six different absolute magni- 
tude ranges. 

The inverse concentration index, given by C = 
petroR50/petroR90, has been used in our analysis. For typ- 
ical galaxies, C ranges from 0.3 (concentrated) to 0.55; for 
comparison, a uniform disk would have C = 0.75. 

In the distribution of the control galaxies (Controls A; 
Figure [5} , the peak of the red sequence gets slightly red- 
der (0.75 to 0.90) and more concentrated (0.38 to 0.32) with 
brighter magnitude, while the blue sequence also gets redder 
(0.45 to 0.60) but at an approximately constant concentra- 
tion index of 0.45. At fainter magnitudes ( — 18 > M r > —19) 
there are few control galaxies on the red sequence, and there 
is just a blue, low-concentration cloud. At brighter magni- 
tudes, it is the blue sequence that begins to disappear, with 
only the red sequence clearly defined. The AGN population 
peaks in approximately the same place as the peak of the 
red sequence for the control sample, hence also gets redder 
and more centrally concentrated with brighter magnitude. 

Figure [6] shows the colour-concentration relations for 



the sample of composite galaxies, along with Controls B. 
The distribution of Controls B varies in the same way as 
Controls A in Fig. [5] but the composite galaxies show a 
slight difference to the Type 2 AGN. The peak of the 
composite population again becomes more centrally con- 
centrated with increasing magnitude, but there is almost 
no change in colour. A more complex morphology in the 
peak of the composite distribution also emerges in the range 
—20 > M r > —23 [Fig. |6jc-e)], which is inconsistent with a 
simple transition population implied by Fig. [3jb) . 

4.4 Colour— Environment Relations as a function 

Of Mr 

The relationship between galaxy colour and environment 
can also add constraints to the processes involved in galaxy 
evolution. In this section we investigate the variation of the 
projected neighbour density E, with g — r colour and abso- 
lute magnitude. 

Envir onmental densities were determined for SDSS 
galaxies bv lBaldrv et al.l (120061 ). using the projected distance 
to the JVth nearest neighbour that is a member of a density 
defining population (DDP). The DDP were galaxies with 
M r < —20. For each galaxy, E was determined from the 4th 
and 5th nearest DDP galaxies within Azc — 1000 krns -1 . 
The environment measurements were only reliable to z < 
0.085 so our sample sizes were reduced to 22 747, 21 666, 
44 592 and 42 966, for Type 2 AGN, composites, Controls A 
and Controls B, respectively. This also included a modest re- 
jection of galaxies whose density measurements were highly 
uncertain. 

Figure [7] shows how mean galaxy environment varies 
with colour and absolute magnitude. The top row shows the 
AGN sample and Controls A, while the bottom row shows 
the composites and Controls B. We have plotted the mean 
log E, in colour bins of 0.05, for three magnitude ranges. We 



© 2007 RAS, MNRAS 000. 011121 



8 P. B. Westoby et al. 





Figure 5. Colour vs. concentration index for different luminosity ranges. The contours are determined in 0.05 colour X 0.05 concentration 
bins. The contour levels are on a logarithmic scale, starting at 64 and doubling each contour [Plots (a) and (f) had insufficient numbers 
for contours]. The solid red lines/red points are Type 2 AGN, and the dashed blue lines/blue points are 'Controls A' galaxies. 
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Figure 6. Colour vs. concentration index for composite sample, in stellar-mass bins of 1 M r . The solid green lines/green points are 
composites, and the dashed blue lines/blue points are 'Controls B' galaxies. See Fig. [5]for details. 
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have also further split these up into two concentration bins; 
C < 0.40 (nominally bulge-dominated galaxies; solid lines), 
and C > 0.40 (nominally disc galaxies; dashed lines). 

From Fig. |4j a simple colour divide between the red 
and blue sequences occurs at (g — r)p e t ro = 0.7. Below 0.7 
lie mainly late-type galaxies, and above 0.7 are the early- 
type galaxies. Generally speaking then, Fig. [7] shows that 
late-type galaxies occur in less dense environments, while an 
increase in colour across to the red sequence coincides with 
more dense environments. Comparing just the solid lines 
in Fig. [7] — i.e., bulge-dominated systems — suggests that 
on the red sequence, inactive galaxies occur in more dense 
environments than active galaxies — Type 2 AGN or com- 
posites. This effect is also stronger at fainter magnitudes. In 
contrast, galaxies with C > 0.4 show little difference in en- 
vironment density across the different classes of galaxy, but 
all still show an increase in density with increasing colour. 

At bluer colours, the error bars are greater as there 
are generally fewer late-type galaxies. Since AGN are also 
predominantly red, there is little reliable data for this class 
at (g - r) Pet ro < 0.6. 

Figure [8] shows the distributions of logE, for galaxies 
with —19 > M r > —22. The top panel consists of galaxies 



with (g — r)p et ro > 0.7 and C < 0.4, so predominantly 
concentrated red galaxies, while the bottom panel shows the 
distribution for galaxies with 0.5 < (g — r)p etro < 0.7 and 
C > 0.4. There is little or no difference in the distribution 
of log E across the various galaxy classes for the bluer, disc 
galaxies (bottom panel), with all populations favouring less- 
dense environments, with a peak at log E = —0.5. 

For the concentrated red galaxies, however, this is not 
the case. As in Fig. [7] we again see that AGN hosts favour 
less dense environments, but this time we also see a slight 
bimodality in the environments of inactive control galax- 
ies. There is one peak, which falls in line with the peaks of 
the AGN and composite populations, which is again at log 
E = —0.5, but there is also a secondary peak at a higher den- 
sity of approximately log E = 0.75, which most likely corre- 
sponds to a unique population of galaxies in higher density 
environments. 



5 DISCUSSION 

Having conducted a robust classification of 'active' and 'in- 
active' galaxies based on their emission line properties and 
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log S 

Figure 8. Top panel: Normalised distribution of log E for galaxies 
with (g — r)p etro > 0.7 and C < 0.4 (red sequence). Bottom 
panel: Normalised distribution of log £ for galaxies with 0.5 < 
(9~ r )petro < 0.7 and C > 0.4 (blue sequence). Solid red lines are 
Type 2 AGN, and solid blue lines are Controls A, while dashed 
green lines are composites, and dashed blue lines are Controls B. 

locations in the BPT diagram, we have created an AGN cat- 
alogue, with a carefully selected control catalogue; and a cat- 
alogue of composite galaxies with another separate sample 
of well-matched control galaxies. The control galaxy samples 
both display a bimodal distribution in rest-frame colour, and 
also Ha line strength. There is an AGN sequence that peaks 
along the red sequence and extends down to the red edge of 
the blue sequence. Composite galaxies appear to be slightly 
less luminous, and bluer, with the peak of the distribution 
possibly in the valley betw een the red and blue sequences 
(see also lMiller et alJl200&t . 

There are a number of possible explanations for this 
result: one interpretation is that star formation maintains 
nuclear activity — i.e., starburst win ds are fuelling the 
central black hole ijNandra et al.ll2007t ). However, the fact 
that the majority of our AGN lie on the red sequence, 
suggests a shorter timescale for star formation associated 
with nuclear activity. This then points in the direction of 
AG N activity suppressing star formation in m assive galax- 
ies l|Springel et al.ll2005l ; Hopkins et"aH 120051 ). In this sce- 
nario, there needs to be an initial injection of gas to feed 
the AGN, and then over time the AGN actually terminates 
star formation. In the local universe, AGN feedback has been 



suggested as a mechanism to e xplain the observed effects of 
star formation quenching (e.g.. iBower et al]|2006l ). Gas ex- 
pelled during the feedback is believed to cool and relax the 
central regions of the galaxy resulting in passive evolution 
of the host galaxy, with stellar winds fuelling the AGN. This 
theory could also account for the observed M. — a relation. 

This would then naively predict that galaxies in the 
process of AGN feedback would lie in the valley between 
the blue and red sequences, as they gradually reduce the 
amount of on-going star formation until they end up on the 
red sequence. However, the efficiency of AGN feedback must 
be considered here, as it is unlikely that every AGN is able 
to suppress star formation in the galactic bulge every time 
it ignites. It may be that feedback can only occur during a 
quasar phase, and not during periods of weak AGN activ- 
ity. If this is the case, then the weak AGN would not be 
responsible for the galaxy appearing in a 'transition' region, 
and the star-formation, AGN and dynamical timescales in 
the nucleus are likely to be the important factors for these 
galaxies. 

While we do find many galaxies in this 'valley', particu- 
larly in the composite sample, it is difficult to conclude that 
this is a transitioning population, because of the complex 
morphology of composites in the colour-concentration plane 
shown in Fig. [6] Composites appear to be either already on 
the red sequence, or on the red edge of the blue sequence 
with a distribution in concentration index that is consistent 
with them remaining in the blue population. Other than 
stellar populations informing us whether or not there is a 
black hole, there is no obvious correlation between colour- 
con centration and th e abili ty to fuel a black hole. 

iDrorv fc Fisherl l|2007l ) explain galaxy bimodality in 
terms of two different types of bulges — classical bulges, 
which are dynamically hot systems, and "pseudobulges" , 
which are dynamically cold, disc-like structures. The dif- 
ferent structures are then thought to be indicators of the 
galaxy's evolution. Classical bulges form through violent re- 
laxation during major mergers in early epochs, when the 
environment density was high, making them red today. In 
contrast, pseudobulges are considered disc-only galaxies that 
have not undergone a major merger since forming their disc, 
hence are much younger, and bluer. Fig. [5] supports this 
interpretation, and shows that AGN are generally classi- 
cal bulge-dominated systems (E to Sa). This result, along 
with the M. — a relation indicates that the SMBH correlates 
with the classical bulge component, and not the galaxy as 
a whole, although the role of the host galaxy in triggering 
and fuelling the nuclear activity can not be ignored. 

Whilst Fig. HHshows that the majority of composites are 
also bulge-dominated systems, there is also evidence for an- 
other population as well, which is again inconsistent with 
a transition population. However, making a typical blue- 
sequence, disk-dominated galaxy red, does not make it a 
'red-sequence' galaxy, as the latter are typically more mas- 
sive, more concentrated, and of earlier Hubble-type, so it 
cannot be assumed that AGN activity in a less massive, 
late-type galaxy, is turning it into a red-sequence galaxy. 

It has long since been established that early-type (red- 
sequence) galaxies are more cluste red than late- type (blue- 
sequence) galaxies (jPresslerl [l980), but in addition to this 
result, we now find that active galaxies, despite being as red 
as inactive early-type galaxies, do not occur in such dense 
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environments as their inactive counterparts. In other words, 
the fraction of AGN amongst the red-concentrated galaxies 
decreas es in denser environm ents. This result appears in con- 
trast to iMiller et ahl (|2003h who found that the AGN frac- 
tion was constant with environment. However, the combi- 
nation of the passive-galaxy AGN fraction decreasing while 
the passive fraction increases could maintain the constant 
AGN fraction over all galaxies. We also note that we use a 
signif icantly l arger sample with carefully matched controls. 

iLi et al.l (|2006l ) studied the clustering of narrow-line 
AGN in the SDSS, and on scales between 100 kpc and 1 Mpc, 
conclude that AGN are preferentially located at the centres 
of dark matter halos. If this is the case, then a proportion 
of the inactive red-sequence galaxies could be explained as 
satellite galaxies in high density environments, which gives 
rise to the implied increased density shown in Fig. [7] and 
the secondary peak in Fig. |HJa) . Fuelling of AGN could be 
restricted in high-density environments except for galaxies 
at the centres of dark matter halos. 



6 CONCLUSIONS 

We have carried out a robust classification of 'active' and 
'inactive' galaxies in the SDSS, based on their emission line 
properties and locations in the BPT diagram. We have com- 
pared the active galaxies with well selected control galaxies 
matched in redshift, absolute magnitude, aspect ratio and 
radius. We find that: 

• Type 2 AGN host galaxies occur mainly on the red se- 
quence of the CMD, but show increased levels of Ha flux in 
their emission compared to inactive red-sequence galaxies 
(Fig. Ha) and Fig. Ha)); 

• A separate class of composite galaxies appears to peak 
on the blue edge of the red sequence on the CMD, whereas 
the peak of the Ha distribution places composite galaxies 
firmly in the valley between the blue and red sequences 
(Fig. EDO and Fig. Hb)); 

• Colour-concentration relations, however, show a more 
complex, possibly double, morphology in the peak of the 
composite distribution, rather than being in a valley (Fig. [6]); 

• AGN (and composites) are found in less dense environ- 
ments on average then matched inactive red-sequence galax- 
ies. The more clustered inactive galaxies are likely to be 
satellite galaxies in high-density environments (Fig. [7] and 
Fig. |U) . 

The key to understanding this in more detail now lies 
in the dynamics of the central regions of galaxies, to under- 
stand what activates some galaxies, but not others. 



7 FURTHER WORK 

The motivation for the construction of the parent sample 
presented here is an ongoing detailed study of the 3D distri- 
bution and kinematics of gas and stars in active and inactive 
galaxies using the IMACS-IFU on the Magellan telescope. 
Equal numbers of low-redshift (z < 0.05), high- luminosity 
Type 1 Seyferts and Type 2 Seyferts (which includes com- 
posites) are carefully selected by inspection of their SDSS 



spectra to ensure that they have strong emission lines, par- 
ticularly in [OIII] . Control galaxies are then matched to the 
AGN in a similar way to the method presented in §3.2, but 
instead of two controls being produced, 10 to 15 potential 
matches are produced, and the best is selected following fur- 
ther inspection the SDSS images and spectra of the controls. 
This work is an extension of the SA URON-IFU s t udy o f 
nearby active and inactive galaxies by [Dumas et all (2007). 

The study of the characteristics of a large statistical 
sample (this paper) complements detailed kinematic studies 
of a well-selected, but necessarily smaller sample of galaxies, 
which can identify the physical mechanisms that drive the 
trends observed in the statistical sample, through detailed, 
spatially resolved studies of individual galaxy nuclei. 
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